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1 Abstract 

^ : 

t-H | The purpose of this paper is to investigate the fundamental problem of the non-uniform 

lO ■ subsonic motion of a point force and line forces in an unbounded, homogeneous, isotropic 

| medium in analogy to the electromagnetic Lienard-Wiechert potentials. The exact closed- 

form solutions of the displacement and elastic fields produced by the point force and line 
forces are calculated. The displacement fields can be identified with the elastodynamic 
Lienard-Wiechert tensor potentials. For a non-uniformly moving point force, we decompose 

' the elastic fields into a radiation part and a non-radiation part. We show that the solution 

of a non-uniformly moving point force is the generalization of the Stokes solution towards 
the non-uniform motion. For line forces the mathematical solutions are given in the form of 

■ time-integrals and, therefore, their motion depends on the history. 
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1 Introduction 



An important item in elastodynamics is concerned with the radiation and the waves produced 
by the non-uniform motion of body forces. This is a fascinating and interdisciplinary research 
topic. The radiation problem has attracted the interest of researchers from different fields 
such as applied mathematics, material science, continuum mechanics, and seismology (see, e.g, 
[1, 19, 12, 2, 20]). A fundamental question is: what is the elastic radiation caused by non- 
uniformly moving point forces? 

In elastostatics, the so-called Kelvin problem concerns with the displacement and elastic 
fields produced by a static point force. In elastodynamics the displacement field generated 
by a time-dependent concentrated point load was first presented by Stokes [22] (see, e.g., [1, 
2, 11]). In the Stokes problem, the body force is considered as a concentrated load of time- 
dependent magnitude. The Stokes solution can be considered as the first mathematical model of 
an earthquake [4]. Concentrated line forces with time-dependent magnitude were studied by de 
Hoop [6] and Achenbach [1]. The wave-motion caused by a line force moving non-uniformly 
in a fixed direction was considered by Freund [10]. A non-uniformly moving line force in an 
anisotropic elastic solid was studied by Wu [23]. 

The radiation problem of point forces is three-dimensional so that Huygens' principle prevails. 
Using the Helmholtz decomposition, the so-called retarded potentials were given for the waves 
produced by body forces in elastodynamics (see, e.g., [1, 19]). A more general expression for the 
retarded potential in elastodynamics was given by Hudson [12]. Elastodynamic fields propagate 
with finite velocities. There always is a time-delay before a change in elastodynamic conditions 
initiated at a point of space can produce an effect at any other point of space. This time-delay 
is called elastodynamic retardation. 

In electrodynamics, radiation is caused by the non-uniform motion of an electric point charge. 
The electric and magnetic potentials of such a non-uniformly moving point charge are called the 
Lienard-Wiechert potentials. The corresponding electric and magnetic field strengths consist 
of velocity-depending fields and acceleration-depending fields. The last ones are the fields of 
radiation. This is a standard topic in electromagnetic field theory and is covered in a lot of books 
on electrodynamics (e.g. [15, 13]). It is quite surprising that nothing has been investigated in 
this direction in the elastodynamics of moving point forces. No solution of a non-uniformly 
moving point force analogous to the Lienard-Wiechert potential can be found in standard books 
on elastic waves (e.g. [1, 19, 9, 7, 12, 2, 20]). 

The purpose of the present paper is to investigate the fundamental problem of the non- 
uniform motion of a point force as well as line forces in an unbounded, homogeneous, isotropic 
medium in analogy to the electromagnetic Lienard-Wiechert potentials. We consider the sub- 
sonic motion. The paper is organized as follows. In Section 2, we present the framework of elas- 
todynamics and we formulate the equation of motion. In Section 3, using the three-dimensional 
elastodynamic Green tensor, we calculate the elastodynamic Lienard-Wiechert potential of a 
point force. In Section 4, using the Lienard-Wiechert potential of a point force, we determine 
the elastic distortion and the velocity fields (particle velocity) of the medium caused by the 
non-uniformly moving point force. In addition, we specify the radiation fields proportional to 
the acceleration of the point force. The limit to the Stokes solution is performed in Section 5. 
The static limit of the displacement and elastic fields of the non-uniformly moving point force 
is given in Section 6. In Section 7, using the two-dimensional Green tensors, we give the gen- 
eral solution of the two-dimensional non-uniformly moving line forces. We close the paper with 
conclusions in Section 8. 
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2 The elastodynamic equation of motion 

In elastodynamics [11], the force balance law reads 1 

Pi - djOij = Fi , (1) 

where p, a and F are the linear momentum vector, the force stress tensor and the body force 
vector. In the theory of linear elasticity, the momentum vector p and the stress tensor a can be 
expressed in terms of the physical state quantities, namely, the velocity vector (particle velocity) 
v = u and the elastic distortion tensor f3 = (grad u) T of the medium which can be derived from 
a displacement vector u by means of the following constitutive relations 

Pi = pvi=pui, (2) 
&ij = Cijki Pki = Cijki d[U k , (3) 

where p denotes the mass density and Cijki is the tensor of elastic moduli. The tensor C^i 
possesses the following symmetry properties 

Cijki = Cjiki = Cijik = Ckhj ■ (4) 

If we substitute the constitutive relations (2) and (3) in Eq. (1), we obtain the force balance law 
expressed in terms of the displacement vector u 

pui - CijkidjdiUk = Ft. (5) 

The solution of Eq. (5) can be represented as a convolution-integral in space and time. In 
an unbounded medium and under the assumption of zero initial conditions, which means that 
u(r,to) and u(r,to) are zero for to — > — oo, the solution of u reads 



/t f'OO 
/ Gijir-r^t-t^Fjir'^dr'dt'. (6) 
-oo J —oo 

Here, Gij is the elastodynamic Green tensor of the anisotropic Navier equation defined by 

[Sik pdtt - Cijkidjdi] G km = S im S(t)S(r) , (7) 

where 5(.) denotes the Dirac delta function and 5ij is the Kronecker delta. The tensor of elastic 
moduli for isotropic materials is given by 

Cijki = ^5ij$ki + p{5ik8ji + SaSjk) , (8) 

where A and fi are the Lame constants. Substituting Eq. (8) in Eqs. (5) and (7), we obtain 
respectively the isotropic Navier equations for the displacement vector 

[Sij p d tt - Sij n A - (A + p) didj] uj = Fi , (9) 

and for the elastodynamic Green tensor 

[Sij pd tt ~ Sij pA - (A + p) didj]G jm = S im S{t)S{r) , (10) 



1 Spatial differentiation is denoted by dj = d/dxj, and for the differentiation with respect to time t we use the 
notation pi = dtPi- 
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where A denotes the Laplacian. 

When the material is isotropic and infinitely extended, the three-dimensional elastodynamic 
Green tensor reads [18, 9, 2, 20] 

G " {r ' t) = ^r P \7% i(t ~ r/CT) + ? (4 5{t ~ r/CL) " k 5{t ~ r/CT) ) 

+ (^-%)i£%*(t-r)d,}, (11) 

where r = yxf + x\ + x| . It should be pointed out that the tensor in Eq. (11) is the retarded 
Green tensor. Here cl and ct denote the velocities of the longitudinal and transversal elastic 
waves (sometimes called P- and S-waves). The two sound- velocities can be given in terms of the 
Lame constants (ct < cl) 



cl = \ , ct = a /-. (12) 

V P V P 

The elastodynamic Green tensor (11) is a tensor with support along the two sound-cones r = opi 
and r = c\,t as well as in between them. Moreover, it should be mentioned that it consists of 
near-field and far-field terms. The first two terms in Eq. (11) decay as 1/r and, thus, they 
are the far-field terms. The last term in Eq. (11) decays more rapidly as 1/r 2 which gives the 
near- field term (see, e.g., [5]). 

3 The elastodynamic Lienard-Wiechert potential of a point force 

Now we consider the non-uniform motion of a point force of total strength Qj(t), situated at 
the position s(t). Then the point force is 

Fi(r,t) = Qi(t)S{r - s(t)) for r £ I 3 , teR. (13) 

Moreover, only subsonic source-speeds will be admitted (| V\ < ct). Substitution of Eq. (13) in 
Eq. (6) and integration in r' lead to 

Ui{r,t)= f G l3 {r-s{t'),t-t')Q 3 {t')dt' . (14) 

J — oo 

The structure of the Green tensor (11) produces in Eq. (14) three characteristic integrals which 
we have to calculate. If we substitute the elastodynamical Green tensor (11) in Eq. (14) and 
use the relation 

i rr/cT /-1/ct 

/ T5(t-T)dr= / n5(t-Kr)dK, (15) 

r Jr/c L A/cl 

where k is a dummy variable with the dimension of slowness = 1/ [velocity], we obtain 
1 /■* f Sij 5(t - t' - R{t')/c T ) 



c\ R(t>) 

+ (i. S (t - t> - R(t')/c L ) -^S(t-t'- R(t')/cr)) 



(16) 
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where R(t') = r - s(t') and R(t') = [R m (t')R m (t')] 1/2 = \r - s(t')\. Now the time-integration 
in Eq. (16) can be performed. We express the integrals in terms of retarded variables by appeal 
to the relation [21, 3] 



/ 



S(f(t'))g(t')dt' = 



9(f) 



|d//df| 



at/(t')=0 



(17) 



Mathematically, the factor l/|d//dt'| is the Jacobian of the transformation from t' to the new 
integration variable /(<')■ This mapping between the two variables is one-to-one if the Jacobian 
is different from zero. A sufficient condition for this is that the velocity of the source (point 
force) is less than the slowest wave speed (|V| < ct). The first integral of Eq. (16) can be carried 
out with 



/ 



5(t-t'-R(t')/cr)Qj(t') , 
R(t') 



Qj(t') 



R(t>) - V m (t')R m (t')/c T 



(18) 



t'=t T 



where we introduced the so-called transversal retarded time = t'(r,t) which is the solution 
of the condition 



t-t' - R(t')/c T = 0. 



(19) 



V = s is the velocity of the moving body force and R is the distance vector from the position 
of the point force s, the sender of elastic waves, to the point of the observer r, the receiver of 
the elastic waves. The second integral is 



/ 



Ri{f)Rj(f)Qj(f) 
R?(t') 



5(t -t'- R(t')/c L:T ) dt 1 



Ri(t')Rj(t') 



Qj(t') 



R 2 (t>) R(t<) - V m (t')R m (t')/c L , T 



t'=*L,T 

(20) 



Here £l = tf(r,t) denotes the longitudinal retarded time which is the solution of the equation 



t-t' - R(t')/c L = 0. 
We perform the third integral as follows 



(21) 



/( 



3^(0^(0 _ A_)q^) J^nS(t - t> - «fl(0)d«df 



R 3 (t') 

r 1 /^ fm t (t')R 3 (t') 

hi 



/i/ CL V R 2 (t<) 

with t K = t'(r,t) as solution of the equation 

t — t' — k R(t' 



Qj(t') ndn 



R(f) - KV m (t')Rm(t') 



0. 



t'=t K 



(22) 



(23) 



The retarded time t K is an effective retarded time for the ^-integration with limits (1/cl, V c t)- 
The solutions of t^, t^, t K are unique if |V| is less than ct- Thus, for subsonic motion the 
solutions of Eqs. (19), (21) and (23) are unique. The retarded times are a result of the finite 
speeds of propagation for elastodynamic waves. In Eqs. (18), (20) and (22) we have used the 
relation 



dt' 



t'=t Iet 



V m (t')R m (t') 
cR(t>) 



t tr, 



>0 for|F|<c T , c = c t ,c l ,1/k and i ret = t T ,th,t K , 

(24) 
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where f(t') = t - if - R{t')/c. 

Thus, carrying out the integration in t' in Eq. (16), we find the explicit expression for 
the displacement field of a non-uniformly moving point force which we call the elastodynamic 
Lienard-Wiechert potential of a point force 



Ui(r,t) 



1 1 



4irp 1 c T 



+ 



13 R\t) 



QAt') 



R(t>) - V m (t')R m (t')/c T 
Qj(t') 



t'=t T 



LL Ri{t>) R{t')-V m (t')R rn (t>)/c h 



rl/cT 

+ J Ann 



R 2 {t') 



f=*L 



R(t')-KV m (t')R m (t')\ 



(25) 



t'=t K 



where R(t'), V(t') and Q{t') are to be evaluated at the corresponding retarded times. The first 
Jjj-term in Eq. (25) has the form of the acoustic Lienard-Wiechert potential [3, 8] and if Q is 
time-independent, this term reduces to the form of the well-known electric Lienard-Wiechert 
potential of a point charge in electrodynamics (see, e.g., [15, 13]). Due to the appearance of two 
velocities of the elastic waves, the elastodynamic Lienard-Wiechert potential of a point force has 
a more complicated but rather straightforward structure. Eq. (25) consists of three characteristic 
pieces. The first term is the transversal one, transmitting with speed ct, and it corresponds 
to ,5- wave motion. The second term is the longitudinal one, transmitting with speed cl, and 
it corresponds to P-wave motion. The third term is neither longitudinal nor transversal and it 
gives contribution arriving at speeds between the two characteristic ones. This shows that this 
factor represents a combination of P-wave and 5-wave motion. 

It is also important to note that the following elastodynamic Doppler factors appear in 
Eq. (25): 

1 - ni(t rc t)Vi(t TCt )/c for |V| < c T , c = c T ,c L , l//c and t ret = tT,t^,t K , (26) 
where rij = R4/R. 



4 Radiation of the elastic fields of a point force 

In order to derive the elastic fields, we have to calculate the gradient and the time-derivative of 
the displacement field (25). The appearance of the retarded times (19), (21) and (23) in Eq. (25) 
makes the differentiation more difficult. To evaluate these derivatives we need some more basic 
derivatives (see also [3]). We introduce 



P c (t') = R(t') -V m (t')R m (t')/c, 
and obtain after a straightforward calculation 

dtn \R k {t') 



c = ct, cl, 1/k 



dx k 



t'=t ie t 



[cPc(t')\ 



(27) 



(28) 



t'=t ret 



d k [Q 3 {t')] 



t'=«ret 



dt' dQjjt' 
dxk dt' 



cP c (t>) 



Qj(t') 



(29) 



t'=iret 
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dk 



[Pc(t')\ 



*'=tret 



V m {t')R m {t') 



R k (t') + 1 



F 2 (0 



R k (t') - 



Pc(t') 



V k (t') 



t'=*ret 

(30) 



and 



<9fc 



rp^OiMOl 




wo / 


P 2 (0 


*'=tret 


U 2 (0 V 



^ + fM] | ^(0 ( 5 h + MM) 



c p c (0 y i? 2 (to 

2R i (t')R j (t')R k (t') 



cPc(f) / 



R 3 (t')P c (t') 



(31) 



t'=t r . 



where R(t'), R(t'), V(t'), P c (0 and Q(0 are to be evaluated at the corresponding retarded 
times. Using Eqs. (29)-(31), we find for the gradient of the displacement field (25) 



Sij 



fli(0^ (0^(0^(0 



R 2 {t') 



p|(0 



1 

+ ? 

t'=t T C L 



^(0^(0^(0 Qj(f) 



R 2 {t') 



pm 



+ 



+ 



+ 



Q,(0 



rl/c-r 
+ / dKK 2 
■/l/CL 

Ri(t')Rj(t') 



?,R t {t')R 3 {t') \R k (t')Q j (t') 



P 2 (0 



V p * 2 (0 



t'=t« 



12,(0 



R 2 (t')P T (t') 



R 2 {f) ypKov 



Pi(0 



orPr(tf) / RHt'W) 



^•(OWA 2R i (t')R j {t')R k {t l ) 



c T Pr(0 



R?(t>)P 2 (t>) 



t'=t T 



Qj(t') 



R 2 {t') Pl(t') 



R k (t') +(l- 



CL 



p,(i' 



vs(0**(0 



R 2 {t')PL{t>) 

rl/c T 

+ / dnnQjit') 
Ji 



Ri(t') 



c L P L (0 ) R?(t')P h {t')\ 5jk+ c L P L (0 



Vj(i?)R k (if)\ , 2R i (t')R j (t')R k (t') 



+ 



'l/CL 



3Ri(t')Rj(t') 



\ 1 / 



<5ifc + 



+ 



3Pj(0 
P 2 (0^(0 
6fli(0flj(0fl]fc(0 

P 3 (0^(0 



^(0^(0 

p 2 (0 % J (^(0iU0^(0 

+ (i - K 2 F 2 (0)P fc (0 - «p«(0^(0 

/cVK0flfc(0\ 3^(0 / /e^(f)l? fc (0 
^— • 1 - -~ - ■■ - - A ^'fc H ^— 



*'=t L 



p K (f ) ; R 2 (t>)p K (t>) v 



p k (0 



(32) 



This is the elastic distortion field produced by a non-uniformly moving point force. For a 
constant strength Qj = constant, the ^j-term in the i-p-expression agrees with the corresponding 
one given in [3] for the gradient of the Lienard-Wiechert potential in acoustics. The elastic 
distortion consists of parts depending on the velocity V, parts depending on the acceleration V 
(elastic radiation part), and parts depending on Q. Note that the dot over V and Q indicates 
a derivative with respect to the argument, namely the corresponding retarded time. 

In order to calculate the time derivative of the displacement (25), we also need the relations 



~8i/_ 
~dt 



t'=*ret 



-Ritr 

PS'). 



(33) 



t'=tret 



7 



9t[Qj(t')]\ t , 



t'=t r 



dt' dQjjt') 
~dt 



ill' 



t'=t r , 



R(t' 



Pc(t 



t'=t Tet 



(34) 



dt 



[Pc(t')\ 



t'=u, 



(V m (t')R m (t f ) - V 2 (t')) + V m (t')R m (t') 



(35) 



t'=t iet 



and 







R 2 {t>) 


*'=tret 



[R(t f )Pc(t 



^(^(0^(0 + ^(0^(0) 
2^(0^(0^(0^(0 

R 3 (f)P c (f) 



(36) 



*'=W 



Using Eqs. (34)-(36), we obtain for the time derivative of the displacement field (25) 



Vi(r,t) 



-M- 

4np | c\ 



SijR(t') - 



^(0^,(0 WO 



R(f) 



P|(0 
SijRit'] 



rl/cT 

+ I Ann 

Jl/CL 



SRi^Rjit') 



+ 



1 



t'=t T L L 



R t (t')R 3 (t>) Qjjf) 
R(f) Pl(t')\ 



*'=t L 



+ 



Qj(0 



+ 



Q;(0 



% - ^^y ) ([^(0^(0 - ^ 2 (0] ^ + MOiUO 
1 (v<4i\R a'\ > T/^^?? ^ 2^(f)i? j (t / )y m (t / )i? m (t / ) 

+ i?(o^(o ( i( } j( } j( w )} «o 

" ^o^(o ( i( } j( } j( } l( )} — — 



t'=t T 



t'=t h 



1/cT 



1/CL 



dnnQj(t') 



3i?i(0-Rj(0 

^7) 



- ^)^7y([Kn(0^m(0 - ^ 2 (0]«#(0 + Hn(0#m(0) 



fl(0^(* 



3 -(^(f )Ji,-(0 + ^(0^(0) + ^ W(OK>(0^(0 



WW) 



t'=t« 



(37) 



This is the velocity field (particle velocity) produced by a non-uniformly moving point force. 
For a constant strength Qj = constant, the <5y-term in the ii-expression agrees with the corre- 
sponding one given in [3] for the time derivative of the Lienard-Wiechert potential in acoustics. 
Again the velocity vector consists of parts depending on the velocity V, fields depending on the 
acceleration V (radiation part), and fields depending on Q. The dot over V and Q indicates 
again a derivative with respect to the argument, namely the corresponding retarded time. 



5 The Stokes solution as limit of a non-uniformly moving point 
force 

In the current section, we show that the Stokes solution is contained in our general results for 
a non-uniformly moving point force. In this sense, our solutions (25) and (32) are the correct 
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generalizations of the Stokes solution towards the non-uniform motion. If the position of the 
point force is fixed, which means that s is time-independent and therefore V = 0, we recover 
from the displacement (25) the famous Stokes solution of a concentrated point force with time- 
dependent magnitude (e.g. [11, 9]) 



Ui(r,t) 



1 



AnpR 
+ 



'IJ 



3RiRj 



"1/ct 
1/cl 



\ rl/c T 

$ij I / nQj(t — kR) dn 



(38) 



The first terms in Eq. (38) are usually called the far-field terms since they behave as 1/R and 
the last term in Eq. (38) is called the near-field term (see, e.g, [2, 20]). From Eq. (32) and 
after some mathematical manipulations, we find the corresponding displacement gradient of the 
Stokes solution (e.g. [11, 9]) 



1 



4np 



bRiRjRk SijRk + SjkRi + SikRj 



i? 3 



\ rl/c T 

I / nQj(t — kR) dn 
/ Ji/cl 



+ 



+ 



+ 



QRjRjRk _ SjjRk + SjkRj + SjkRj 
SijRk 



R 3 



1 



Qj(t-R/cL) 



j- Qj(t - R/ct) 



i? 3 



R 



Qj(t - R/ct) + — Qj(t - R/c T ) 

CT 



RiRjRk 
R 4 



\ Qj(t - R/cl) - -i- Qj(t - R/ct) 



(39) 



Thus, Eqs. (25) and (32) give the correct Stokes solution as limit. It can be seen that the 
Q-terms in Eq. (39) behave as 1/R (far-field terms or radiation terms) and the Q-terms behave 
as 1/R 2 (near- field terms). 



6 Static limit of a non-uniformly moving point force 



In this section, we give the static limit of the Lienard-Wiechert potential and the elastic distortion 
of a point force as a further check of the obtained results. 

For the static limit, we set V = and Q = constant and substitute Eq. (12) and A = 
2^i//(l — 2v) in Eq. (25). If we perform the integration in k and arrange in proper order the 
appearing terms, we recover the displacement field of the Kelvin problem (see, e.g., [11]) 



Ui(r) 



Qj 



16tt/x(1 - v)R 



(40) 



where R = r — r' and v is Poisson's ratio. From Eq. (32), we find in the static limit the 
displacement gradient of the Kelvin problem (see, e.g., [11]) 



Afc(r) = - 



Qi 



167T/i(l - v)R 3 



(3 - 4v)5ijRk - SikRj - SjkRi + 



3R{RjRk 
R~ 2 



(41) 



Thus, in the static limit, we recovered the displacement and displacement gradient fields of 
the Kelvin problem, which is concerned with a concentrated point force in three-dimensional 
elastostatics. 
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7 The elastodynamic Lienard-Wiechert potentials and elastic 
fields of non-uniformly moving line forces 

We proceed to derive the two-dimensional Lienard-Wiechert potentials and the elastic fields 
of non-uniformly moving line forces. We consider line forces moving non-uniformly and with 
time-dependent magnitude. The line forces are parallel to the X3-direction. In two dimensions, 
the in-plane line force is given by 

F a = Q a (t)5(R(t)), (42) 

and the anti-plane line force reads 

F 3 = Q 3 (t)6(R(t)), (43) 

where R(t) = r — s(t) £ M 2 and a = 1,2. The field variables are independent of the x 3 - 
component. 

If the material is infinitely extended and isotropic, the two-dimensional elastodynamic Green 
tensor of plane-strain reads [9, 14] 

G^r, *) = f H { t - r /cL) - H<t - r/or) 



- ¥ (^-ryctH(t - r/cL) - ' H(t - r/c T )) 

\ Jt 2 -r 2 /c^ ' 

and the elastodynamic Green tensor of anti-plane strain is given by 

1 H(t- r/or) 
27T P C T ^t 2 - r 2 /c\ ' 



(44) 



, ^ 1 H(t-r/c T ) , . 

G 33 (r,t) = - T } 1 T> , (45) 



where H(.) denotes the Heaviside step function and r = ^x 2 + x\. 

If we substitute Eqs. (44) and (42) in Eq. (14) and perform the integration in r', we find for 
the displacement field 



/ \ 1 
u a (r,t) 



2vrp 



L w ] { #m sjF) + v #m ) Sh{t V dt 

. ( R a (t')Rp(t>) ( R a (t>)Rp(t')-6 a pR 2 (t>) \ t 2 \ . 

L W ] { R\t>) S ^ } + { rW) ) SrM ) dt . 



(46) 

The notation here is 

t = t-t', S 2 (t') = t 2 -^1, Sl{t') = t 2 -^P-. (47) 

c T c L 

The two retarded times ir and th are the roots of S^(t') = and S^(t') = 0, respectively, which 
are less than t. The solving of the conditions S^(t') = and S 2 (t') = is non-trivial for a 
general motion and can be complicated. For a subsonic motion, the solutions for the retarded 
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times £t and are unique. The corresponding elastic distortion and velocity of a non-uniformly 
moving line force are 



— !— 8y 

2vrp 7 



/»'•■(/'; ' V|iM - i fifj I— 



"J S T (f ) 



(48) 



and 

v a (r,i) 



/■tL 

■/ — oo 

W) 

-oo 



R a {t')Rp{t>) I- f //„(/')// ,■(/') tf„,/*Vi } 



R A {t') S L (t') 



R A {t') 



Ra(t')Rp(t') ^ + ( R a (t')R[](tJ-6 a0 R 2 (t>) ^ g [[/ 



i? 4 (t') 



i? 4 (t') 



V s T (f 



(49) 



Now we consider the anti-plane line force. If we substitute Eqs. (45) and (43) in Eq. (14) 
and perform the integration in r', we find for the displacement field of a line load of body forces 
pointing in the X3-direction 



u 3 (r,i) 



1 



t T 



s T (t>) 



dt 1 . 



The corresponding elastic distortion and velocity fields are 



-oo 



Q3(t') 

5 T (f) 



dt'. 



v 3 (r,t) 



*T 



Q 3 (0 



dt'. 



(50) 



(51) 



The displacement fields of non-uniformly moving line forces play the role of the two-dimensional 
elastodynamical Lienard-Wiechert potentials. The two-dimensional Lienard-Wiechert poten- 
tials (46) and (50) are time-integrals over the history of the motion and, thus, they are charac- 
terized by an afterglow. Such two-dimensional wave motion possesses a 'tail' characteristic for 
the so-called diffusion of waves. For that reason line forces are haunted by their past. Neverthe- 
less, for the two-dimensional Lienard-Wiechert potentials (46) and (50) we are left to evaluate 
time-integrals of considerable complexity, which only in some simple cases yield results of ele- 
mentary functions in a closed form. Also, the calculation of the retarded times is not a trivial 
task. The mathematical complexity of the integrals (46)— (51) is the same as of the integral 
expressions of non-uniformly moving straight dislocations given by Lardner [16] and Lazar [17]. 

If we put s(t') = in Eqs. (46) and (50), we recover the displacement fields of concentrated 
line forces given by de Hoop [6] and Achenbach [1]. For the in-plane line force, de Hoop [6] and 
Achenbach [1] used a more complicated but equivalent representation of the Green tensor (44). 



8 Conclusion 

Exact analytical solutions of the displacement and of elastic fields have been calculated for 
point and line forces moving non-uniformly in an unbounded, elastic, isotropic body. We have 
investigated the subsonic motion (|V| < or)- We have shown that the displacements can be 
interpreted as the elastodynamic Lienard-Wiechert potentials caused by body forces. For a 
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point force, we calculated explicitly the radiation parts of the elastic fields. We have proven 
that our solution of a non-uniformly moving point force is the correct generalization of the 
Stokes solution and of the solution of the Kelvin problem. In the case of line forces the Lienard- 
Wiechert potentials are given in the form of time-integral representations, which cannot be 
further simplified for the general non-uniform motion. 
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